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Abstract. We discuss mesoscopic effects in quantum dots, nanoparticles and nuclei. In quantum dots, we focus on the 

statistical regime of dots whose single-electron dynamics are chaotic. Random matrix theory methods, developed to explain the 

Cu I statistics of neutron resonances in compound nuclei, are useful in describing the mesoscopic fluctuations of the conductance 

i-^ , in such dots. However, correlation effects beyond the charging energy are important in almost-isolated dots. In particular, 

jy^ ■ exchange and residual interactions are necessary to obtain a quantitative description of the mesoscopic fluctuations. Pairing 

^ ' correlations are important in metallic nanoparticles and nuclei. Nanoparticles smaller than --^ 3 nm and nuclei are close to the 

fluctuation-dominated regime in which the Bardeen-Cooper-Schrieffer theory is not valid. Despite the large fluctuations, we 

find signatures of pairing correlations in the heat capacity of nuclei. These signatures depend on the particle-number parity of 

protons and neutrons. 

B 

-6 : 1- INTRODUCTION 

O I Mesoscopic systems can include a large number of particles yet are sufficiently small that finite-size quantum effects 

are important. A good example of a mesoscopic system is a quantum dot, a sub-micron scale conducting device in 

which up to several thousand electrons are confined by electrostatic potentials 1 1]. Another example is a nanoparticle 

- an ultrasmaU metallic grain |21. The discreteness of the energy spectrum of a quantum dot or a nanoparticle becomes 

^ ' important at the low temperatures achieved in laboratory experiments. 

The transport properties (i.e., conductance) of quantum dots and nanoparticles can be measured by connecting them 
to leads. Linear conductance measurements provide information on the ground state of the system, while properties of 
■«^ , excited states can be extracted from non-linear conductance experiments. Furthermore, by changing a gate voltage, it 

f^ ' is possible to control the number of electrons in the device. 

^^ I To a certain extent a compound nucleus can also be considered as a mesoscopic system. Several methods and ideas 

^P ' introduced in nuclear theory have indeed inspired developments in mesoscopic physics. 

In this article we review some of the mesoscopic effects in quantum dots, metallic nanoparticles and nuclei. We 
consider some analogies and differences between quantum dots and nuclei and between nanoparticles and nuclei. 
In small vertical dots, the confining potential is harmonic, leading to shell structure and the phenomenon of magic 
Fj^ . numbers familiar from atomic and nuclear physics [IJ. The spectra of such dots can usually be understood in terms of 

r^ ' shell effects plus an exchange interaction. 

Q , For larger lateral dots, the symmetries responsible for the shell structure are broken and the single-electron dynamics 

O ■ in the dot become mostly chaotic Olj] . In this regime, the main interest is in understanding the statistical properties of 

the spectrum and conductance sampled from dots with different shapes or magnetic fields rather than in explaining the 
precise properties of a particular dot. In this regime random matrix theory (RMT) 1 5, 6 j has proven to be very useful. 
. ^ RMT was developed by Wigner, Dyson, Mehta and others in the 1950's and 1960's to explain the statistical properties 

S i of the neutron resonances in compound nuclei. 

The simplest model of a quantum dot is the constant interaction (CI) model in which the electron-electron interaction 
is taken to be the classical charging energy. In this limit the conductance is simply related to the single-particle wave 
functions in the dot, and the RMT approach is justified by the chaotic dynamics associated with the single-particle 
Hamiltonian. However, in an almost-isolated dot, residual interactions beyond the charging energy are important and 
must be taken into account to achieve a quantitative understanding |3, 4|. Furthermore, the spin degrees of freedom 
beyond the simple degeneracy factor become crucial. The exchange interaction associated with the spin degrees of 
freedom has explained some of the discrepancies between the experiments and the predictions of the CI model. 
Understanding the interplay between mean-field behavior and residual interactions is at the frontier of research in 
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mesoscopic physics, making this field even more attractive to many-body theorists. 

Pairing correlations play an important role in metallic nanoparticles |2|l- Bulk metals become superconducting at 
low temperatures, a phenomenon explained by the Bardeen-Cooper-Schrieffer (BCS) theory fT\. BCS, as a mean field 
theory, is valid in the limit where the single-particle mean level spacing is small compared with the pairing gap. 
However, in nanoparticles smaller than ^ 3 nm, this condition no longer holds and the grain is in a crossover between 
the BCS limit and the regime dominated by fluctuations of the order parameter The nucleus also belongs to this 
crossover regime since the single-nucleon mean level spacing within a shell is typically smaller than the pairing gap. 
The effects of pairing correlations are much more subtle in this crossover regime |2] and it is apriori not clear whether 
any thermal signatures of the bulk pairing transition survive in this regime. We expect to see some similarities between 
nanoparticles and nuclei in the manifestation of pairing correlations. 

We remark that in comparison with natural systems such as nuclei, quantum dots have the advantage that their 
parameters (e.g., size, shape, number of particles, coupling to the outside) can all be easily controlled by the experi- 
mentalist. This makes quantum dots a practical tool to study finite-size effects in small quantum systems. 

Quantum dots are discussed in Sec.|2] where we focus on the statistical regime of almost-isolated ballistic dots. The 
CI model is briefly discussed in Sec. 12. II The application of RMT to the mesoscopic fluctuations of the conductance in 
the CI model is described in Sec. 12. 21 The CI plus RMT model can explain only some of the data, and its limitations are 
summarized in Sec. 12. 31 Sections lT?l through l2.6l discuss correlation effects beyond the CI model. In Sec. l2.4l we review 
the so-called universal Hamiltonian, obtained in the limit of a large number of electrons. In addition to the charging 
energy term, it contains an exchange interaction that depends on the total spin S of the dot. In Sec. 12. 51 we present a 
closed expression for the conductance in the presence of such an exchange interaction. This is a non-trivial many-body 
problem but through the use of spin projection it can be reduced to quantities that characterize non-interacting spinless 
electrons. The remaining residual interaction terms are small. However, they play an important role as is explained in 
Sec.Ol 

Nanoparticles are briefly discussed in Sec. |3 emphasizing the significance of pairing correlations in finite-size 
systems. Section^discusses pairing correlations in hot nuclei. In Sec. l4.1l we explain how thermal properties of nuclei 
can be calculated beyond the mean field approximation using quantum Monte Carlo methods. In Sec. l4.2l we use these 
methods to study signatures of pairing correlations in the heat capacity of finite nuclei. We find that these signatures 
depend on the parity of the proton and neutron numbers. 



2. QUANTUM DOTS 

An example of a quantum dot fe] is shown in Fig.^ A quantum dot is typically made by forming a two-dimensional 
electron gas at the interface of a seminconductor heterostructure (e.g., GaAs/AlGaAs). Applying electrostatic poten- 
tials to metallic gates the electrons are confined to a 2D sub-micron region yj. 
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FIGURE 1. Quantum dot used in the experiment of Folk et al. FS|. Right: a scanning electron micrograph of the dot (top view). 
Left: a schematic drawing. 

The transport properties (i.e., conductance) of a quantum dot are measured by connecting it to leads. At sufficiently 
low temperature (below ^100 mK), the coherence length L^ of the electrons in the dot is larger than the size of the 



dot L and the system is called mesoscopic. Phase decoherence of the electron results from inelastic scattering from 
other electrons or phonons. 

The electron in the dot scatters elastically from impurities with a mean free path /. The dot is diffusive when / ^ L. 
The dot is ballistic when it has little disorder, i.e., Z > L. In the latter case, transport is dominated by scattering from 
the boundaries of the dot. 

In small vertical dots (up to A^ ^ 20 electrons) the confining potential is often harmonic and leads to shell structure 
19]. The addition spectrum (i.e., the energy requires to add an electron into the dot) shown in Fig.|2]displays maxima at 
magic numbers A^ = 2,6, 12 that correspond to fermions with spin 1/2 in a 2D harmonic potential. Secondary maxima 
are observed at half-filled shells A^ = 4,9, 16. This is a result of Hund's rule from atomic physics - the electrons are 
added with parallel spin up to mid-shell so as to maximize the (attractive) exchange interaction. 
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FIGURE 2. Shell structure observed in the addition energy of a small vertical dot. The dot has a shape of a disk with a harmonic- 
like confining potential. Maxima are observed for filled shells of a 2D spin-degenerate harmonic oscillator (magic numbers) and 
secondary maxima correspond to half-filled shells (following Hund's rule). From Ref. j^l. 

In larger dots the symmetries responsible for shell structure are broken and the single-particle dynamics are mostly 
chaotic (i.e., unstable with respect to the initial conditions) 1 10]. Figure|3]shows a schematic drawing of a ballistic dot 
with a chaotic trajectory scattered from the boundaries of the dot. This regime is the statistical regime, in which we are 
interested in understanding the mesoscopic fluctuations of the conductance |3]. The relevant time scale in a ballistic 
dot is the time Tf to cross the dot. The statistical regime has some analogies with the compound nucleus, in which 
equilibration in (n, 7) reactions is achieved through collisions of the incoming neutron with other nucleons (known as 
the "doorway" state mechanism). Here the relevant time scale is the equilibration time Teq. 




FIGURE 3. A schematic drawing of a ballistic quantum dot attached to two leads. A typical trajectory scattered several times 
from the boundaries of the dot before exiting. 

The width of a resonance in a quantum dot is controlled by the coupling between the dot and the leads. Denoting 
by Tc the transmission coefficient for the electron to decay into channel c in a lead and by Th the recurrence time (i.e., 
the Heisenberg time), the resonance width is given by 
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where we have used T// = h/d (d is the mean level spacing in the dot). 



An open dot is usually characterized by a large number of open channels with T^ ^l soF^ d,a situation analogous 
to the compound nucleus in the regime of overlapping resonances. The conductance of an open dot exhibits Ericson 
fluctuations flT] versus a gate voltage. 

As we pinch off the point contacts, effective tunneling barriers are formed in the dot-lead interfaces and 7^ <C 1 . 
In such almost-isolated dots, F <C A and the conductance exhibits a series of peaks as a function of gate voltage. An 
almost-isolated dot has some analogies with the regime of isolated neutron resonances in the compound nucleus. The 
top panels of Fig.|6lbelow compare the conductance peaks in a quantum dot with the neutron resonances observed in 
the total cross section to scatter thermal neutrons from heavy nuclei. 

2.1. Constant interaction (CI) model 

The distance between the conductance peaks in a quantum dot is typically much larger than the single-particle mean 
level spacing d. This is explained by a simple model, known as the constant interaction (CI) model. In the CI model, 
the Coulomb interaction energy is taken to be the classical electrostatic energy e^N^ /2C, where A^ is the number 
of electrons in the dot and C is the capacitance of the dot. The single-particle levels in the dot are spin-degenerate 
(ff = ±1) and depend only on the orbital label A. The total Hamiltonian of the isolated dot in the CI model is 

V-, ^ t e^N^ 

Ha = l^{£x - eaV^)a\^axa + -^ ■ (2) 

The term linear in the electron-number operator N = Lact '^l^'^'Acj describes the effect of the gate voltage Vg (a is a 
constant determined by the ratio between the gate-dot capacitance and the total capacitance of the dot plus gate). 

2.1.1. Coulomb blockade 

In general the tunneling of an additional electron into the dot is suppressed by the Coulomb repulsion of the 
electrons already in the dot, a phenomenon known as Coulomb blockade. By changing the gate voltage of the dot, 
it is possible to compensate for this Coulomb repulsion. The dot will conduct at its degeneracy points, i.e. when 
Ef + Ef,s{N — 1) — Egfi{N). Here gp is the Fermi energy of the electron in the leads and EgaiN) is the ground-state 
energy of a dot with A^ electrons (we have assumed the temperature to be sufficiently low so that only the ground state 
of the dot participates in the conductance). Using Eq. (|3, we find that the effective Fermi energy gp = Ep +eaVg at 
the degeneracy points satisfies 

e,= {N-\yl + e,. (3) 

The conductance displays a series of peaks at values of the gate voltage given by (|3j- These are known as Coulomb- 
blockade peaks. The spacing A2 between two successive peaks is given by 

A2 = Agp = eVC + (e^+i - En) . (4) 

The charging energy e^/C is typically much larger than the mean level spacing d and therefore A2 ~ e^/C is 
approximately constant. The phenomenon of Coulomb blockade is demonstrated in Fig.0] 



2.1.2. Rate equations approach 

When a typical width F of a resonance in the dot much smaller than both the mean level spacing d and the 
temperature (i.e. F ^ d,kT), the coherence between the dot and the leads can be ignored. It is then possible to use 
a rate equation approach to calculate the conductance through the dot flZ|. For low temperatures kT <C d, only one 
resonance level A (near the Fermi energy) contributes to the conductance. The conductance peak is broadened by 
temperature and the important information is coded in the peak height G. This peak height is given by G 0= ^ wSX 
where 

_ 2 r(^r^ 
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eaVg = (N+l/2)^ 
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FIGURE 4. A schematic drawing of Coulomb blockade in a quantum dot: (a) when eaVg = Ne /C there is a charging energy 
gap on both sides of the Fermi energy that prevent the tunneling of an electron into the dot; (b) when eaVg = {N + l/2)e^/C there 
is no gap at the Fermi energy and an additional electron will tunnel into the dot. 



Here F, is the partial width of the level X to decay into the left (right) lead. 



2.2. Conductance fluctuations and random matrix theory 



For a chaotic dot, the conductance peak height exhibits large fluctuations as a function of the gate voltage (see, e.g., 
in Fig.|6jc)). In Ref. 1 13] it was suggested that these fluctuations can be described quantitatively using random matrix 
theory (RMT). RMT was introduced by Wigner 1 14] in the 1950's to explain the fluctuation properties of the neutron 
resonances in the compound nucleus, and further developed by Dyson 1151 . Mehta and others. The Hamiltonian matrix 
is assumed to be random except that it must satisfy the underlying space time symmetries. 

In RMT we have an ensemble of Hamiltonians H described by a Gaussian distribution law 
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(6) 



where a^ is a typical variance of a matrix element of H. There are three types of ensembles: (i) Gaussian orthogonal 
ensemble (GOE) of real symmetric matrices, (ii) Gaussian unitary ensemble (GUE) of complex hermitean matrices, 
and (iii) Gaussian symplectic ensemble (GSE) of real quaternions. These three ensembles are labelled by j3 = 1,2 
and 4, respectively. The GOE corresponds to conserved time-reversal symmetry while the GUE corresponds to broken 
time-reversal symmetry (e.g., in the presence of an orbital magnetic field). The GSE is relevant in the presence of 
spin-orbit scattering and will not be discussed here. 

In nuclei, RMT is used to describe the statistical fluctuations of the many-body levels and eigenstates at finite 
excitation energies Ill6iri7ll . This approach was justified by the complexity of the nuclear Hamiltonian. In chaotic 
quantum dots, RMT is used to describe the statistical fluctuations of the single-particle levels and eigenstates. This is 
based on the Bohigas-Giannoni-Schmit conjecture ] 18] stating that the statistical quantal fluctuations of a classically 
chaotic system are described by RMT. The universality of RMT holds for time scales f > T, where 



T/ (chaotic ballistic dot) 
(compound nucleus) . 



(7) 
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The ballistic Thouless energy of a quantum dot is the energy scale associated with T/, i.e., Et = f'/'^f- The 
dimensionless Thouless conductance gj is then defined as the number of single-particle levels within a Thouless 
energy window 

Ej 



Estimating T/ ~ L/vp (where Vf is the Fermi velocity) we find gj 
vaUd within gj levels around the Fermi energy. 



(8) 



'N. According to Q, the RMT description is 



The distributions of the conductance peak heights in almost-isolated chaotic quantum dots were predicted using 
RMT |,13;]. The partial width amplitudes correspond to different components of an RMT eigenvector. Thus each of the 
partial widths P'' is described by a Porter-Thomas distribution of the corresponding symmetry j3 . Assuming the leads 
to be uncorrected, the following closed-form expressions are found for the peak height distributions I.13> ,19.] : 
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where Kq and Ki are modified Bessel functions. The distributions (|9} and ([Toj are universal and depend only on the 
symmetry class. They were measured in two independent experiments |8, 20] and found to agree with the theoretical 
predictions. The comparison between theory and experiment is shown in Fig.|5] 
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FIGURE 5. Conductance peak height distributions in Coulomb blockade quantum dots. Left: the measured distributions of 
Ref. 1 20] for small quantum dots in the absence [panel (a)] and in the presence [panel (b)] of a time-reversal symmetry-breaking 
magnetic field B. Right: the measured distributions of Ref. |8] for large quantum dots. The insets show the same distribut ions 
(histograms) on a log-linear scale. The solid lines in all panels are the theoretical predictions of Ref. 11.311 given by Eqs. J9j and <10> 
in the absence and presence of a magnetic field, respectively 

It is interesting to compare the statistics of the conductance peak heights in an almost-isolated quantum dot with the 
statistics of the neutron resonance widths in a compound nucleus. The top panels of Fig. |6l compare a typical series 
of neutron resonances in a compound nucleus (left) with a typical series of conductance peak heights in a quantum 
dot (right). The bottom panels of Fig. |6l compare the distributions of the neutron resonance widths (left) with the 
distribution of the conductance peak heights in the absence of magnetic field (right). 

Other signatures of quantum chaos in the conductance peak height statistics which we do not discuss here are the 
parametric peak height correlations (e.g., at different values of the magnetic field) 12211 and the increase of the average 
conductance with orbital magnetic field 12311 . They were measured in Refs. [8] and [24^, respectively. 



2.3. Limitations of the constant interaction model 



While the CI plus RMT model was successful in explaining (at least qualitatively) some of the observed mesoscopic 
fluctuations of the conductance there have also been major discrepancies. One such discrepancy was observed in the 
statistics of the peak spacing A2 ]25, 26, 27, 28]. At low temperatures, when only the ground state contributes to the 
linear conductance, the peak spacing is given by 



A2 =Eg,{N+ I) +Eg,{N- I) -2Eg,{N) , 
where Eg^iN) is the ground-state energy of the dot with A^ electrons. 
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FIGURE 6. Neutron resonances in a compound nucleus [panels (a) and (b)] and conductance peak heights in a chaotic quantum 
dot [panels (c) and (d)]: (a) total cross section of n+ Th as a function of the incoming neutron energy E„ (in eV); (b) the 
distribution of the normalized neutron resonance width f = F/F on a log-linear scale. The histogram is the data of 223 resonances 

12111 and the solid line is the Porter-Thomas distribution /'(f) <x f ^'' e^ ' ; (c) a series of Coulomb blockade peak in a quantum 
dot versus gate voltage Vg at zero magnetic field 1;^; (d) the distribution of the normalized conductance peak heights G = G/G. The 
histogram is the data of 600 peak heights and the solid line is the RMT prediction of Ref. I13ll . From Ref. y|l- 



In the CI model, A2 ~ e^/C + Ae, where Ae is the spacing between two successive single-particle levels [see Eq. 
0]. The charging energy e^/C is a constant and the statistics of A2 are determined by Ae. Since the single-particle 
levels are spin-degenerate, we expect the peak spacing distribution to be bimodal, i.e., a superposition of a 5-function 
(corresponding to Ae = 0) and the Wigner nearest-neighbor spacing distribution (derived from RMT). Including noise, 
we expect a bimodal distribution shown by a solid line in the inset of Fig. [7] However, the observed distribution of 
Ref. 1I27II shown in Fig.0 is not bimodal and its shape is closer to a Gaussian. The absence of bimodality is a common 



observation of all the corresponding experiments I25L 126112 71 12 




FIGURE 7. The experimental peak spacing distribution -P(Vg) where v^ = (A2 — A2)/A2 fTf\. The histogram is the data 
containing 10800 peaks at S f^ and kT ~ 0.8d, and the solid line is a Gaussian fit. The inset is the peak spacing distribution 
in the CI plus RMT model before (dotted line) and after (solid line) convolution with a Gaussian noise. From Ref. 12711 

The peak spacing statistics has also been studied as a function of temperature |27]. The observed width (7(A2) of 
the spacing fluctuations (symbols in Fig.llO> is seen to be suppressed when compared with the CI plus RMT model 



(long dashes). More detailed studies of the peak height statistics at finite temperature |29l] also reveal significant 
discrepancies. The kT = Q.\d experimental distribution, shown by the shaded histogram in the left panel of Fig. E] 
suggests that the CI plus RMT model (dashed histogram) overestimates the probability to find small peak heights. 
Furthermore, the ratio C7(^niax)/gmax between the width and average of the peak height distribution is suppressed in 
comparison with the predictions of the CI plus RMT model (see right panel of Fig. II 11 1. 



2.4. The universal Hamiltonian 



The origin of the discrepancies discussed in Sec. l2.3l is in the oversimplified treatment of electron-electron interac- 
tions by the CI model. In an almost-isolated dot, charge is quantized and it is necessary to consider correlations beyond 
the simple charging energy term. Of particular importance are correlations that are associated with the spin degrees of 
freedom. 

In the absence of interactions, the energy is minimized in a state of lowest spin, i.e., S — Q for an even number of 
electrons and 5=1/2 for an odd number of electrons. Any configuration with a higher spin requires the promotion of 
an electron to a higher single-particle level and therefore an increase in energy. Such a state of minimal spin (shown 
on the left in Fig.|S} is known as a Pauli state and is typical of metals. 
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FIGURE 8. Left: Pauli state, typical of metals. Right: Hund state, typical of atoms. 

The Coulomb interaction is repulsive and is minimized in a spatially antisymmetric state. Such a state requires a 
spin-symmetric state to ensure antisymmetrization of the total wave function, i.e., a state of maximal spin. This effect 
can be described by an effective exchange interaction which is minimized in a state of maximal spin. In atoms, such 
a state is known as the Hund state (shown on the right in Fig.|S). In general, the ground-state spin of a quantum dot 
is determined by the competition between the total non-interacting energy of the electrons in the confining potential 
and the exchange interaction. In chaotic quantum dots in which shell structure is absent, the ground state is often 
intermediate between a Pauli state and a Hund state. 

In a chaotic dot, the randomness of the single-particle wave functions induces randomness in the interaction matrix 
elements. In general one can separate the electron-electron interaction into an average and fluctuating parts. The 
standard deviation of the fluctuating part behaves as ^ d/ gj 1.3(1 .31.1 for a diffusive dot and ^ d^/lngr/gr for a 
ballistic dot I32l 133113411 . In the limit of a large Thouless conductance (i.e., for a large number of electrons) only the 
gT'-independentpart of the average interaction survives. This leads to the so-called universal Hamiltonian 1351 141 

i/universal = J^is^ ~ eaVg)al^axa + ^ - JsS^ + dpiJ.T^T , (12) 

where S = jY^allaa' ^aa^aa'^aa' is the total spin of the dot (a is the vector of the three Pauli matrices) and 
T' = Y.ia'^ aj^_ creates pairs of spin up-spin down electrons. Js = v'a/3;6a is the average strength of an exchange 
interaction matrix element while Jc = v„„.nft describes an average Cooper channel interaction. The Cooper channel 
term is allowed only in the orthogonal symmetry (i.e., in the absence of an external magnetic field). In quantum dots, 
it is repulsive and does not lead to any instability |4|]. In the foUowing discussions on quantum dots we will ignore this 



term. However, in nanoparticles and nuclei the pairing interaction is attractive and plays an important role, as we shall 
discuss in Sec.|3]andE] Note that the universal Hamiltonian reduces to the CI model when J^ = and Jc = 0. 

The statistical fluctuations of the single-particle wave functions and energies within a Thouless energy band are still 
described by RMT. The new important term in the universal Hamiltonian is the constant exchange interaction, and the 
various statistical properties of the dot will depend on the particular value of the exchange interaction constant J^. RPA 
estimates provides an expression for /( as a function of the gas constant r^. A typical value of r^ for the dots used in 
the experiments is rj ~ 1.2 and the corresponding RPA estimate is 7j ^ 0.3d. 

In the universal Hamiltonian plus RMT model, there is a certain probability for the quantum dot to have spin S in 
its ground state L36..34.1 . In Fig.|9]we show these ground-state probabilities for different spin values as a function 
of Js (measured in units of d). In the CI plus RMT model the ground state for an even number of electrons is 
always characterized by 5 = 0. However, for a realistic value of the exchange interaction J^ = 0.3 d, the ground-state 
probability of spin 5 = 1 is about 20%. 
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FIGURE 9. Probability of the ground state to have spin 5 versus the exchange constant Js (measured in units of d). Solid (dashed) 
lines correspond to the orthogonal (unitary) symmetry. From Ref. I34ll . 



2.5. Conductance in the presence of interactions 

The conductance formula (|5} is based on a non-interacting model (except for the charging energy term which 
is constant for a fixed number of electrons). In the presence of interactions, it is necessary to generalize the rate 
equation approach using the many-body wave eigenfunctions of the dot. This was accomplished in Ref. |37]. The 
input parameters are the energies of the many -body eigenstates of the dots with A^ and A^+ 1, respectively, and the 
transition widths Y^'- between the A^-electron dot in an eigenstate / and the A^+ 1 -electron dot in an eigenstate j 
involving an electron tunneling from the left (right) lead. In the linear response regime it is generally necessary to 
solve a linear set of equations but in special cases a closed formula can be derived L37il 
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(13) 



Here /),- is a Fermi-Dirac distribution for an energy-conserving transition between the many-body states / and /' and 
P\ is the equilibrium probability of finding the dot in the eigenstate / of A^ electrons. 



2.5.1. Conductance formula for a constant exchange interaction 



In particular, a closed formula of the form ( I13> holds when the orbital occupations ni = ni^ + nx^ are good 
quantum numbers 137, 331 . This is true for the universal Hamiltonian \\2\ since [S^,n;L] ~ 0. However, Eq. \\2>\ 
requires a summation over all the many-body eigenstates of the dot with A^ and A^ + 1 electrons. The number of 



these eigenstates increases combinatorially with the number of single-particle levels, and ( I13> becomes impractical at 
higher temperatures. Using spin projection techniques we were able to perform the combinatorial sum over all orbital 
occupation numbers in a closed form I38ll . 

The eigenstates of J12> (with Jc = 0) are characterized by the orbital occupation numbers n — (ni ,n2i ■ ■ ■) of the 
single-particle levels. The manifold of all Slater determinants with m given singly occupied levels is of dimension 
2"'. In the absence of exchange, all these 2'" states are degenerate with energy Y.X^X^'^X ("ot including charging 
energy). States with good spin S can then be constructed by taking the appropriate linear combinations of these 
2'" states. The problem is equivalent to the coupling of m spin- 1/2 particles to total spin S. In general, there are 
several degenerate states with the same spin S and they can be distinguished by additional quantum numbers 7. The 
eigenstates of good spin S and spin projection S, = M are thus labeled by \N'a'YSM) . Their corresponding energies are 

Applying spin (and particle-number) projection, the linear dimensionless conductance can be written as 

g = £(wf +w('))g,, (14) 

X 

where gx are the single-level conductances (JSJ. The thermal weights w\_ and w^ collect the contributions from 
processes in which an electron is added to an empty level or to a singly occupied level, respectively. They are given by 

>^f =4^/7;l,^,,P^.s Y. (25'+l)/(e^,5), (15) 

s y=s±i/2 



w 



X 



^Y.'^x.n+x.s'Pn+x.s' E (25+l)[l-/(es%)], (16) 



y 5=S'±l/2 

L _AN+\) AN) ~ ™ 
'S-^n'S' ^^n5 ^ ^F' ^^^ qi 

given by 



where g^^ = £^,„, — £^^^ — gp. The quantity Pj^^s is the probability to find the dot with A^ electrons and spin S. It is 



P^5 = e-/^fe+f/~.sl/Z, (17) 

where Fj^s = — jS^'lnTr^r.^e^"^^"'^^"^!' is the free energy of A^ non-interacting electrons with total spin 5, Un.s = 
e^N^ /2C — JsS{S + 1) — EfiN and Z is the grand-canonical partition function restricted to A^ and A^+ 1 electrons. The 
quantities b^f^^ — j{{hi — l){nx—2))i^s ^nd c^^^^ = j{ni{nx — l))N,s ^6 projected quantities at fixed particle- 
number A^ and spin S. 

The free energy Fn,s and the coefficients b,c can be expressed in closed form in terms of the free energy F^ and 
canonical occupations {nx)g of q spinless fermions with energies e;^: 



g-pfN-s^g ' Y+^ T-^-e T+*+' T-(-^+'). (18) 

CX.N.S — _R/F., J_F., ^ _ftfF,, iP., \ • ^^"> 



^A w 5 is found is by replacing fii in Eq. ( I19> with {I —fix)- 

For a finite number of single-particle levels, the quantities Fq and {ni)y can be expressed in closed form using 
particle-number projection [see, e.g., Eqs. (140) in Ref. p|]. 

2.5.2. Peak spacing statistics 

In Ref. p81, we have used the conductance formula of the universal Hamiltonian (discussed in Sec. 12.5. H to 
calculate the statistics of peak spacings and peak heights. For a chaotic dot, the single-particle spectrum and the partial 
widths Fx are calculated from random matrices with the appropriate symmetry. For each sample, the conductance is 
calculated as a function of Ef using the formulas of Sec. 12.5.11 and then maximized to find the peak height. FigureFTOI 
shows the width C7(A2) of the peak spacing fluctuations as a function ofkT/d. The solid line corresponds to Js — 0.3 c/ 
and is in excellent agreement with the experimental results (symbols). We note that this value of /s is not fitted but is 
extracted from the value of the gas constant r^ of the experimental samples using an RPA estimate. When compared 
with the CI model {J^ = 0, long dashed), we observe that the exchange interaction suppresses the fluctuations of the 
peak spacing. 
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FIGURE 10. Standard deviation of peak spacing fluctuations ct(A2) in the presence of magnetic field as a function of kT /d. 
Lines are theoretical results for Js = 0,0.3c/ and Q.5d and the symbols are the experimental results of Ref. 12711 . From Ref. 13^ . 



2.5.3. Peak height statistics 

The peak height statistics I3a.l39ll are shown in Fig. ^2 The left panel shows the peak height distribution at the lowest 
experimental temperature of kT = 0.1 d. The calculated distribution for Js = 0.3d (solid histogram) describes well the 
experimental distribution (shaded histogram). On the other hand, the CI model (dashed histogram) overestimates 
the probability of small peak heights. The ratio C7(gniax)/^max between the average and width of the peak height 
distribution is shown versus kT /d in the right panel of Fig.fH] This ratio is found to be suppressed by the exchange 
interaction. The Js = 0.3d results (solid line) are in good agreement with the experiment (symbols) for low and 
intermediate temperatures kT < 0.6 d. At higher temperatures there are deviations that are still not understood. 
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FIGURE 11. Peak height statistics in chaotic dots (in the presence of magnetic field). Left: the experimental peak height 
distribution |29] for kT = O.lrf (shaded histogram) is compared with the theoretical distributions for 7^ = (dashed histogram) 
and for Js = O.Bd (solid histogram). Right: the ratio (T(gmax)/|max between the standard deviation and the average value of the 
peak height as a function of kT /d. Theoretical results are shown for Js = (long-dashed line), 7, = Q.'id (solid line) and 7s = 0.5 rf 
(short-dashed line). The experimental results of Ref. 1 29] are shown by symbols. From Ref. 
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2.6. Beyond the universal Hamiltonian: residual interaction effects 



For a finite Thouless conductance gr, it is necessary to take into account the residual interactions beyond the 
universal Hamiltonian. In general, the effective Hamiltonian of a chaotic dot in a band of ~ gr levels around the Fermi 
energy is 

L4- i T— 1 -I- -f- 

eafl„^flafj + T 2j Vap-ysaafjOi^^iayaiagf^ , (20) 

aa ^ apyS 



where v is the two-body screened Coulomb interaction in 2D. 

The randomness of the single-particle wave functions induces a randomness of the interaction matrix elements when 
the latter are calculated in the single-particle eigenstates. The possible induced ensembles of random interactions, 
which depend on the underlying space-time symmetries and on features of the two-body interactions, were classified 
in Ref. |40]. The cumulants of the interaction matrix elements are determined by requiring invariance under a change 
of the single-particle basis. The average interaction is characterized by three (two) invariants in the orthogonal (unitary) 
symmetry VaB^yS = ^'o^aySgg +Js5asSHY + JcSaS^8Y^ ^nd has the same form of the universal Hamiltonian (except for 
a linear term in N). The second cumulant for a local two-body interaction in the unitary symmetry is given by 



5v* 



^.ysSvi^v-pa 



-U^{d, 
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l'5/3^5ycj55p) 



(21) 



The corresponding expression in the orthogonal symmetry can be found in Ref. |40]. The constant u^ describes a 
typical variance of a screened Coulomb interaction matrix element, and can be estimated using a contact interaction 
model in the limit of a small screening radius. In a diffusive dot u ^ d/gr 13011 and in a ballistic dot ^ d^lngr/gr 

In the experiments, gr — n^J'N jl varies in the range ~ 30 — 70 and the fluctuations of the interaction matrix 
elements are expected to be small compared with d. The effect of the residual interaction terms can be studied within 
the Hartree-Fock (HF) approximation. Most of these studies has been restricted to spinless electrons f4T,MX?l. 

In a finite dot, screening also leads to a surface charge potential |30] [not included in J20H. The fluctuations of 
the matrix elements of the one-body surface charge potential are of order dj ^fgj, parametrically larger than the 
fluctuations of the two-body interaction matrix elements. However, for the experimental values of gj, the fluctuations 
of the one-body and two-body interaction matrix elements turned out to be comparable. 
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FIGURE 12. Left: the peak spacing distribution at kT — 0.2d for the universal Hamiltonian <12> with J^ — 0.3d (unitary 
symmetry). Right: the peak spacing distribution at kT = 0.22d. The dashed histogram is the calculated distribution for 7, — 0.28c/ 
and includes residual interaction effects in the HF-Koopmans approach. The solid histogram is the experimental distribution of 
Ref. m. From Ref. 01 . 

A prominent experimental feature of the peak spacing distribution is the absence of bimodality |'25, 26, 27, 28]. For 
temperatures below ^ 0.3 d, the peak spacing distribution remains bimodal |33] even when the exchange interaction 
is included (see left panel of Fig. I12t . The absence of bimodality was understood to be an effect of the residual 
interactions f4Z, 43, 33, 34,^4]. Residual interactions (in the presence of spin) were studied using different methods 
including diagonalization for a small number of electrons 1421 . a Strutinsky approach |43], spin density functional 
theory I44ll . and a HF-Koopmans approach 13411 in which the single-particle wave functions are assumed to be 



unchanged upon the addition of electrons | l45l] . In the HF-Koopmans approach, we choose a reference HF many- 
particle state (e.g., an 5 = Pauli state for A^ electrons) and then express the addition energy in terms of the HF 
single-particle levels of the reference state and a few interaction matrix elements. Using the fluctuation properties of 
the HF levels and wave functions we can determine the residual interaction effects on the peak spacing distribution. 
The results (including the effects of surface charge effects and gate voltage scrambling) are shown in the right panel 
of Fig. [2] The theoretical peak spacing distribution compares reasonably well with the experimental distribution. 



3. METALLIC NANOPARTICLES 



The fabrication techniques for quantum dots were extended in the mid 1990's to ultrasmall metallic grains. Figure [T3l 
shows an ultrasmall metallic grain (of size ^^ 2 — 10 nm) connected to leads IMol. For such grains the temperature is 
much smaller than the mean level spacing d and the discreteness of the levels becomes important. By measuring the 
non-linear conductance it is possible to determine the excitation spectrum of an individual grain. 
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FIGURE 13. A schematic drawing of the experimental setup of Ref. l4fll for the transport measurements in Al nanoparticles. 
FromRef. b46i1 . 

In contrast to quantum dots, the metallic nanoparticles become superconducting at low temperatures. Experimental 
and theoretical studies of ultrasmall metallic grains have led to a better understanding of the role of pairing interaction 
in finite-size systems. For a recent review see Ref. |2j]. 

A simple model of a metallic gain is given by 
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(22) 



with J^ < 0. This is essentially the universal Hamiltonian J12I I but with an attractive pairing interaction. We have also 
included a Zeeman energy term in an external magnetic field H (Hb is the Bohr magneton and gs is the ^-factor). 

The pairing interaction in ( I22t has the reduced BCS form. BCS theory explained the phenomenon of superconduc- 
tivity in bulk metals. As a mean field theory it is valid in the limit when the mean level spacing d is much smaller than 
the bulk pairing gap A. In this limit fluctuations can be ignored and a mean field theory is a good approximation. 

This BCS limit holds for metallic nanoparticles that are larger than ~ 5 nm. Indeed a pairing gap was observed 
in the spectrum of such nanoparticles with an even number of electrons. An example is shown in Fig. [3]where the 
excitation spectrum of an Al grain is shown as function of a magnetic field H. A gap is observed for an even number 
(left panel) but not for an odd number of electrons (right panel). In BCS theory, the low excitations for an even number 
of electrons involves the excitation of two quasi-particles, hence a gap of ~ 2 A to the first excited state. On the other 
hand, the lowest states in a grain with an odd number of electrons are described by one quasi-particle states and hence 
no gap is observed between the first excited state and the ground state. 

However, in the smaller grains of size less than ^ 3 nm, we have d > A and the effect seen in Fig.[3]can no longer 
be observed. This is the fluctuation-dominated regime where BCS theory is no longer valid. An interesting question is 
whether signatures of pairing correlations still survive in the crossover from BCS to the fluctuation-dominated regime. 
Theoretical studies have indicated that the crossover is smooth but depends on the particle-number parity (see in |2l 
and references therein). To derive such effects it is important to project on a fixed number of electrons or at the very 
least project on the correct parity of the number of electrons. In Sec. 0]we address pairing correlations in the crossover 
regime by studying properties of nuclei at finite temperature. 
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FIGURE 14. The excitation spectrum of an Al nanoparticle versus a magnetic field H extracted from non-linear conductance 
measurements. Left: transition from odd to even number of electrons. Notice the gap 2A to the first excited state. Right: transition 
from even to odd number of electrons. From Ref. |46]. 

4. NUCLEI 

In a nucleus, a typical spacing d within a shell is smaller than the pairing gap A. Thus the nucleus belongs to a regime 
dominated by fluctuations just as is the case for the smaller nanoparticles. While zero-temperature signatures of the 
pairing interaction are well known in nuclei, much less is known about thermal signatures of pairing in nuclei. 

To find out whether thermal signatures of the pairing transition survive in finite nuclei (despite the large fluctuations), 
it is necessary to take into account correlations beyond the mean field approximation. Interaction effects can be 
described by including fluctuations around the mean field. This is formally described by the Hubbard-Stratonovich 
(HS) transformation |47], in which the Gibbs ensemble e^P^ at inverse temperature j3 is described as a superposition 
of one-body propagators Ua 



tUo 



(23) 



Here Ga is a Gaussian weight and Ua is the imaginary-time propagator describing non-interacting nucleons moving 
in time-dependent external auxiliary fields C7(t). 



4.1. Quantum Monte Carlo methods 

Nuclear properties are well described by the interacting shell model, in which the nucleons move in a mean field 
potential and interact via residual interactions. The model is usually solved by conventional diagonalization of the 
Hamiltonian in the many -particle space. The limitation of this method is the combinatorial increase of the model space 
with the number of nucleons and/or the number of single-particle levels. Shell Monte Carlo (SMMC) methods that are 
based on the HS decomposition were developed to overcome this limitation |48, 49, 50]. The integrand in ( I23> can 
be easily calculated using matrix algebra in the single-particle space. However, the number of auxiliary fields (7(t) is 
very large and the integration in (I23> is done by Monte Carlo methods. 

The thermal expectation value of an observable O is given by 



{0) = 



Tv{Oe-P") _ J^[a]Ga{0)aTvUo 



Trie-P") J^[a]GaTi-Ua 



(24) 



where {0)a = Tr((9t/a)/Trf/(j is the expectation value of the observable O evaluated for a sample a. In SMMC we 
choose the samples according to the probability distribution Wa = Gc7|TrC/(j|, and the thermal expectation value of the 
observable O is then calculated estimated from 



(O) 



La{0)a'^o 



(25) 



where <I>(j = Tr U(j/\Ti Ua \ is the sign of the one-body partition function TiU^- 



Since the number of nucleons is relatively small, it is important to use the canonical ensemble at fixed number 
of protons and neutrons. This can be achieved through particle-number projection 15111 . For example, the partition 
function for A particles in Ns single-particle orbitals is given by 



TtaUo = y e^*'^det l+e*'/''U„ 



(26) 



where ^^ = InmjNs {m = l,...,Ns) are quadrature points, and /i is a chemical potential. 

Often the sign <$>„ fluctuates from sample to sample. When the statistical error of the sign is larger than its average 
value, the statistical errors of observables become too large and the method fails. This is known as the sign problem 
and is generic to all fermionic Monte Carlo methods. A typical effective nuclear interaction suffers from such a sign 
problem. In general the interaction can be decomposed into "good" and "bad" components, depending on whether a 
sign problem is absent or not (for the corresponding component). The dominant collective components of the nuclear 
interaction are good-sign and this property was used in Ref. 14^ to find a practical solution to the sign problem in 
finite nuclei. 



4.2. Pairing correlations at finite temperature 

4.2.1. Hamiltonian 

We have constructed nuclear interactions that contain only the dominant collective components, i.e., pairing and 
multipole interactions. In particular, we have considered the following isospin conserving interaction 1531 

a X 

where 
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and aj,,n.m, = {—Y^'"^^^'"' aj.-m.-m, ( (a similar definition is used forP^^-^O- The mean field single-particle orbitals 
are denoted by a = {nljm) with principal quantum number «, orbital angular momentum I and total angular momentum 
and projection j,m. The corresponding single-particle energies £„ are calculated in a Woods-Saxon potential plus 
spin-orbit interaction |53]. The one-body potential V in (12 8> is the central part of the Woods-Saxon potential, and the 
multipole interaction is obtained by expanding the surface-peaked interaction v(r , r') = —x[dV/dr) {dV /dr')5 (f — f '). 
In the expansion we only keep the multipoles A =2,3,4. The parameter x is determined self-consistently |54| 
X^^ = So" dr r^ {dV /dr) (dp/dr) (p is the nuclear density), and k^ are renormalization factors that take into account 
polarization of the core (we have used k2 —2,kj,^l.5 and ^4 = 1). The pairing coupling strength go is determined 
using experimental odd-even mass differences. All interaction components in ( I27> are attractive and lead to a good- 
sign Hamiltonian. We have studied nuclei in the mass range A ~ 50 — 70 using the Hamiltonian ( I27> in the complete 
.fPS9/2 shell (containing 30 orbitals for protons and 30 orbitals for neutrons). 

The pairing interaction is an important component of the effective nuclear interaction j27l l. In the bulk, BCS theory 
predicts a finite discontinuity in the heat capacity across the pairing transition temperature. However the nucleus is a 
finite system and its heat capacity is expected to behave smoothly as a function of temperature. As discussed above, 
the nucleus is in a regime dominated by fluctuations and an interesting question is whether signatures of the pairing 
transition survive despite the large fluctuations. 

In the absence of multipole interactions, the Hamiltonian Mil has the reduced BCS form, and its eigenvalues can 
be found by solving a closed set of non-linear equations |55]. However, since we are interested in finite temperature 
properties, the SMMC approach offers a more direct approach. Furthermore, the SMMC method allows us to treat the 
more realistic Hamiltonian (|^). 



4.2.2. Heat capacity 

In SMMC, the thermal energy E is calculated as a function of inverse temperature j3 and the heat capacity is 
calculated from a numerical derivative C = dEjdT . However, this leads to large statistical errors in the vicinity of 
the transition temperature (even for a good-sign interaction). We have introduced a novel method to calculate the heat 
capacity that reduces the statistical error by about an order of magnitude Ii561 . This method is based on correlated 
errors. 

The top panels of Fig. [2]show the SMMC heat capacity of ^"^Fe and ^^Fe (symbols) as a function of temperature. 
The calculations are done in the fpgtjn model space using the interaction Ml\ . We observe that the BCS discontinuity 
is strongly suppressed but for the even-mass ^''Fe a 'bump' remains around the neutron pairing transition temperature. 
This bump is correlated with the rapid decrease of the number of 7 = neutron pairs (A^A) (see bottom left panel of 
Fig.ll5>. However, in the nucleus ^^Fe (containing an odd number of neutrons) we do not observe such a bump in the 
heat capacity. 
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Top panels: heat capacities in Fe (left) and in Fe (right). The symbols are the SMMC results in the complete 



FIGURE 15 

fPS9/2 shell and the dotted-dashed lines are obtained in the independent-particle model. For ""Fe we also show the BCS heat 
capacity (solid line). The two finite discontinuities occur at the proton and neutron pairing transition temperatures. Bottom panels: 
the number of nn, pp and pnj = pairs as a function of temperature. The 'bump' in the heat capacity of ^"Fe is correlated with a 
rapid reduction in the number J = neutron pairs. From Ref. 1 56] 

We expect the above signature of pairing correlations in even-even nuclei to be enhanced upon the addition of 
neutron pairs. Figure^]shows the SMMC heat capacities (in the fpg<)/2 model space) for a series of even-mass iron 
isotopes (left panel) and odd-mass iron isotopes (right panel). The bump observed in the heat capacity of the even-mass 
iron isotopes is enhanced as the number of neutrons increases, while for the odd-mass iron isotopes we only observe 
the usual enhancement of the heat capacity with system's size. 



4.2.3. Extending the theory to higher temperatures 



Since only one major shell is included in the SMMC calculations, the heat capacity saturates at higher temperature. 
This happens in the vicinity of the bump (see, e.g., in Fig. I15> . To demonstrate that the bump we observe in even- 
even nuclei is indeed a unique signature of pairing correlations, it is necessary to increase the model space to include 
additional shells. We have extended the shell model theory to higher temperatures by combining the fully correlated 
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FIGURE 16. Systematics of the SMMC heat capacity in iron isotopes. Left: heat capacities of even-mass iron isotopes with 
A = 52 — 62. Notice the enhancement of the 'bump' with increasing number of neutron pairs. Right: heat capacities of odd-mass 
iron isotopes with A = 53 — 61. From Ref. 1561 



partition in the truncated model space with the independent-particle partition in the full space (including all bound 
orbitals and continuum effects) 15711 . 

In the absence of interactions, the grand-canonical many -particle partition function of a nucleus with single-particle 
levels £„/, and scattering phase shifts 5/,(e) is given by 
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(29) 



where 5p{e) = TT^' ^/y(2j + l)d5ij/de is the continuum contribution to the single-particle level density. 

The canonical partition function Z^ at fixed particle-number A can be obtained in the saddle-point approximation 
InZ^ « InZ'-'*-' — PlxA — j In (2;r((AA^)^)) , where ({AA)^) describes the variance of the particle number fluctuation. 

In the presence of interactions we assume 



InZ,, = lnZ,,(r + InZjp — InZ^ 



p.tr ■ 



(30) 



Here Z,, is the partition function in the presence of correlations, Z,,,, is the correlated partition in the truncated 
space, Zsp is the many-particle partition function of the independent-particle model in the full space, and Zsp,tr is 
the independent-particle model partition in the truncated space. In Eq. ( I30t we have subtracted \nZsp,tr to avoid a 
double counting of the truncated degrees of freedom. Eq. ( I30> has the correct form at both low temperatures (where 
InZjp K, XnZsp^tr) and at high temperatures (where correlations become less important). 



4.2.4. Extended heat capacity 



We have used the extended theory of Sec. l4.2.51 to calculate the heat capacity for iron isotopes J57|- The results 
are shown in the right panels of Fig. ^] where the extended heat capacities of ^"Fe and ^^Fe are displayed versus 
temperature (solid squares). In contrast with the truncated SMMC calculations (open squares), the extended heat 
capacity continues to increase monotonically also at higher temperatures. The bump in the heat capacity of the even- 
even nucleus ^"Fe is now clearly observed in comparison with the heat capacity of the independent-particle model 
(dashed line). 

Similar 5-shaped heat capacities were measured in even-even rare-earth nuclei ISal. The experimental heat capacity 
was determined from the many-particle level density p{E) of the corresponding nucleus. This level density was 
measured in heavy ion reactions and the canonical partition function calculated by a Laplace transform Z()3) = 
j dEp{E)e^"^ . The canonical heat capacity was then determined from C = dE /dT where E is the thermal energy 
E = -d\nZ/d^. 
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FIGURE 17. Right: extended heat capacities of Fe (top) and ^^Fe (bottom) versus temperature T J57l]. The extended heat 
capacities (solid squares) are compared with the SMMC heat capacities in the /pgg /2 shell (open squares) and with the heat 

capacities calculated in the independent-particle model (dashed lines). Left: measured heat capacities in '^ Yb and '^'Yb (solid 
lines) are compared with heat capacities in the Fermi gas model (dashed lines) |58]. Even though the experiments are done in a 
different mass region, the shape of the experimental heat capacities are similar to the calculated heat capacities for both even-even 
and even-odd nuclei, respectively. 



5. CONCLUSION 

We have discussed several mesoscopic effects in quantum dots, metallic nanoparticles and nuclei. In particular, we 
have focused on ideas and phenomena in nuclear physics that have found applications in mesoscopic physics and 
nanoscience. 

Random matrix theory, introduced to explain the statistical features of the neutron resonances in compound nuclei, 
turned out to be a useful theory for understanding the mesoscopic properties of chaotic ballistic dots. However, whereas 
in nuclei random matrix theory is used to describe the statistical properties of the excited many -body states, in quantum 
dots it is used to describe the statistical properties of the single-particle spectrum and wave functions. In almost-isolated 
dots, electron-electron interactions beyond the classical charging energy have important effects on the mesoscopic 
fluctuations. A quantitative description can be achieved only when both single-particle chaos and electron-electron 
interactions are taken into account. 

Studies of metallic nanoparticles have shed light on pairing correlations in finite systems. As the size of the 
nanoparticle is made smaller, a crossover occurs from the BCS bulk limit of superconductivity to the fluctuation- 
dominated regime. Of particular interest are pairing correlation effects in this fluctuation-dominated regime. Nuclei 
belong to this crossover regime and can also be used to understand the signatures of pairing correlations in finite 
systems beyond BCS theory. In neutron-rich even-even nuclei we find a 'bump' in the heat capacity that is a signature 
of the pairing transition, while no such signature is found in neighboring nuclei with an odd number of neutrons. 
These theoretical results are in agreement with recent experiments. A similar effect is expected in the heat capacity of 



ultrasmall metallic grains. Theoretical studies in metallic nanoparticles and nuclei suggest that pairing correlations in 
the fluctuation-dominated regime manifest through effects that depend on the particle-number parity. 

Quantum dots, metallic nanoparticles and nuclei are finite-size systems in which the interplay between single- 
particle behavior and many-body correlations leads to a variety of interesting phenomena. We have only described 
some of these phenomena and much remains to be explored in the future. 
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